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Kinetic model of three component, weakly ionized, collisional plasma with a beam of neutral 
^\ , particles is developed. New dispersion relations for linear perturbations are derived and analyzed in 

On ' various limiting cases. 
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I. INTRODUCTION 



It is well-known that neutral beam injection is one of the fundamental fusion plasma heating methods. In general, 

a particle accelerator is used to create fast ion beams (the particle energies are on the order of 100 keV); the ion beam 

^ ■ is then passed through a neutral gas region, where the ions neutralize via charge-exchange reactions with the neutral 

I ' gas. The neutralized beam is then injected into a magnetically confined plasma. Of course, the neutral atoms are 

S I unaffected (not confined) by the magnetic field, but ionize as they penetrate into the plasma. Then the high-energy 

C^ ions transfer fraction of their energy to the plasma particles in repeated collisions, and heat the plasma |jl^-H| . 

^ , In this paper we develop a kinetic model of three component, weakly ionized, collisional plasma with a beam of 

rS i' neutral particles. We employ a kinetic equation for the charged particles of a sort in the weakly ionized plasma with 

^ the Batnagar-Gross-Krook (BGK) model collisional term. Similar model has been developed previously by others 0. 

O ■ In this book authors do not take into account possibility of existence of regular velocity of the neutral particles Q . In 

' ^ I the light of the possible relevance of our model for the heating of plasma by neutral beam injection, we set out with 

J>~j. the aim to generalize results of Ref. M by allowing neutral particles to have regular velocity and seek for possible 

rG ' novelties brought about by this effect. Indeed, the dispersion relations for linear perturbations obtained in this paper 

, ^^ , differ substantially from those of Ref. R . 

In section II we formulate our model and obtain general dispersion relation. In section III we analyze various 
limiting cases of the dispersion relation and discuss the results. 



II. THE MODEL 



> 

o 
o 

We start analysis of the dielectric permittivity (DP) of a coUosional plasma with weakly ionzed, non-degenerate 

On ' plasma when the integral of elastic collisions in the kinetic equation for the charged particles can be apporximated 

On by the BGK term, while it is possible to neglect the collisions between the chaged particles themselvs. Analysis of 

c/3 ' this relatively simple model will be useful for further more complicated case of fully ionized plasma (which, in fact, is 

O , more relevant for the fusion plasma). The latter is beyond the scope of present paper and the separate analysis needs 

C^ ' to be done. 

^~^, The kinetic equation for the charged particles of a sort in the weakly ionized plasma with the BGK model collisional 

1-^ ' term can be written as following [Q: 



Here, Vcm denotes collision frequency of charged particles with the neutrals, which in this model is assumed being 
constant, whereas 



Na= I dp fa, 
and 
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Index a {a = e,i) refers to charged particles (electrons and ions), whereas n — to neutrals. Vq denotes regular, 
uniform velocity of the neutral particles. Finally, Ta is defined by following expression: 
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The specific form of the BGK integral used here is derived from its more general form 



dt 
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where Va/3 is some constant which has meaning of effective collision frequency between particles of a and /3 sort, i.e. 
it characterizes time of momentum relaxation of a sort particles as a result of their collision with particles of /? sort. 
Function $q,^ is determined by following expression: 
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here Vp = {l/Np) J dpvf p. 

It should be emphasized that the BGK coUisonal integral describes accurately collisions only particles of different 
sort. Therefore, it can be used to describe collisions of charged paricles with the neutrals in weakly ionized plasma, 
when the scatteing of charged particles on the neutrals is a dominant process. In the case of fully ionized plasma, in 
spite of its relative simplicity, use of BGK integral is not justified m. 

In what follows, we consider isothermal models of the BGK integral, i.e. we neglect change in temperature of 
charged particles with chage in their corresponding distribution functions. We ought to mention that the results 
obtained here will be qualitatively the same for the non-isothermal model of BGK integral. We further assume that 
the masses and the temperatures of the ions and neutrals do coincide, i.e. rrii = Mn = M and Ti = Tn- In this case 
to the order of '^ m^/M terms we have Te„ = Tg. Thus, in the Eq.(2), under these simplifying assumptions we can 
set Tan — Ta and 



'I'rvn = 
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exp 
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which, in fact, coincides with the Maxwellian distribution function (with the beam having velocity Vq) normalized to 
unity. 

In the static equilibrium state, with the external fields absent, Eq.(l) allows for the only solution /oq — Noa^on- 
In what follows subscript will denote unperturbed and S perturbation of the physical quantities. 

Let us consider small perturbation of the distribution function Sfa which is caused by appearance of small fields E 
and B. After usual linearization of the Eq.(l) we obtain 



dSfa 
dt 
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The solution of the latter equation for the plane monochromatic waves (i.e. E, 5 fa ~ exp 
written as 



-iujt + ik ■ f 
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) can be 
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where rja = {^/Noa) J dpSfa, which is perturbation of the particle number density normalized to equilibrium value of 
the number dinsity. 7]a can be calculated either by integration of the Eq.(7) over momentum or by using the continuity 
equation for the particles of a sort: 
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here, ja denotes charge current of particles of a sort. 

It is known that the complex tensor of DP can be written as 



etj{uj, k) = 5ij -\ o-y (w, k) 

Equj 
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where 5ij is usual Kroneker tensor and crij(aj, k) is the conductivity tensor defined by 

ji = ^iio^ =aij{uj,k)Ej. 



(10) 



In general when e^ tensor is of the type e^ — Sij + AiAj — AiBj, then defining quantities e and e as 
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respectively, we can split tensor from Eq.(9) in the longitudinal and transverse (with respect to wave- vector fc) parts 
as following: 



e,,(c.,fc) = ( % - ^-^ ) £'^(^,fc) + ^£'(c^,fc) 
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Now, inserting expression for (5/q from the Eq.(7) into Eq.(8) and using Eqs.(10)-(12) we obtain following expessions 
for e' and e*'': 
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Here, ujLa = \/{e^Na)/{eQma) and Vto = y/j\Jrn^. The integrals in the Eqs.(14) and (15) may be evaluated by 
choosing the z-axis along k. The integration over Vx and Vy is elementary. Whereas, Vz integral may be expressed 
in terms of a single transcendental function, which called the plasma dispersion function. There are several different 
definitions of this function used in the literature. We use the one given by Melrose |8|: 
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Using Eq.(16) and following intermediate results of integration 
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where, z = {ui + iiyan) / {v^kVTa) , we obtain 

e' = 1 + y ■ ^'" [1 - 0(^)][1 - (fc • K))/('^ + «i^an)] 
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Note, that conventinal kinetic model of three component, weakly ionized, coUisional plasma m is significantly modified 
by taking into account possible existence of a beam of neutral particles. Namely, the expression for the e' is modified 
by additional factor [1 — (fc • Vo)/(w + Wan)]- While the form of the £*'' is not changed by the presence of the beem. 



III. DISCUSSION 



Let us start analysis of the obtained results from longitudinal waves as we have seen that transverse waves do not 
incur any modification by the presence of the beem of neutral particles. The dispersion relation for the longitudinal 
waves reads as following: 
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The latter equation is a transcendental one, thus, in general case, it has many complex solutions Lu{k). Let us consider 
the most interesting ones which correspond to weakly damped oscillations. 

Let us consider, first, high frequency waves, i.e. when uj ^ kvTa, ^an- Using asymptotic expansion for 0(z) [pi 
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we obtain following dispersion relation for the weakly damped waves (Re oj 3> Im lj) 
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(26) 



Here, we neglect the contribution from ions, because it is significant when T^ > Te{M/me)'^ , i.e. when the temperature 
of ions is greater than the temperature of electrons by more than six orders of magnitude. It is unlikely that such 
differences in the temperatures actually do realize in the nature H. Therefore, in the frequency domain concerned, 
the plasma can be considered as a purely electronic, i.e. the role of the ions is reduced only to neutralize the charge 
of electrons. The dispersion relation (26) has to imaginary terms. The first one describes coUisionless Cherenkov 
absorption of the plasma waves. Whereas, the second one has purely coUisional nature and describes dissipation of 
the fields energy in via collisions (electronic friction) H] . The difference induced by the presence of the beam of neutral 
particles is presented by a factor (see, Ref. [Q for comparison) 
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In addition to the high frequency longitudinal oscillations in isotropic coUsionless plasma there also exist low 
frequency oscillations, so called, Ion-acoustic waves. They exist in highly non-isotermal plasma, where Tg 3> T^. 
Phase velocity of these waves lies in the Vn <C oj/k <C Vtc domain. It is obvious, that such waves should also 
exist in coUisional plasma if the collisions are sufficiently rare. Thus, when uj <C h'm and {uj + iven\ ^ kVre in the 
Vti <?C oj /k <^ Vtc phase velocity domain we obtain following dispersion relation 
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In the latter equation we have used also following asymptotic expansion 



(z) = 2z' 
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when 
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to the first order. 

Again, we note that the difference induced by the presence of the beam of neutral particles is presented by a factor 
given by factor Eq.(27) (see, Ref. for comparison). 

This concludes presentaion of the kinetic model of three component, weakly ionized, coUisional plasma with a beam 
of neutral particles. We have generalized the results of Ref. m by allowing neutral particles to have regular velocity 
(i.e. by allowing for the existence of a beam of neutrals). We have shown that the novel, generalized dispersion 
relations for linear perturbations obtained in this paper differ substantially from those of Ref. m . Finally, we would 
like to conclude outlining, once again, the possible relevance of our model for the better understanding of the plasma 
heating process by a neutral beam injection. 
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